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LONG-TERM  GOALS 

The  physical  and  biological  characteristics  of  the  areas  inhabited  by  deep  diving  odontocetes  are 
poorly  understood.  Our  long  term  goals  are:  i)  to  measure  and  characterize  the  biomass  in  areas  and  at 
depths  inhabited  by  beaked  and  spenn  whales;  ii)  to  measure  and  characterize  the  physics  of  these 
environments;  iii)  to  assemble  the  characteristics  measured  (i)  and  (ii)  into  a  depth  integrated,  3- 
dimensional  habitat  model;  the  model  will  include  other  dependent  and  independent  data,  e.g., 
chlorophyll  and  depth,  respectively.  Our  final  long  tenn  goal  is  to  then  apply  the  habitat  model 
produced  to  other  geographic  areas  to  assess  their  likelihood  as  beaked  and  sperm  whale  habitat. 
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OBJECTIVES 


The  past  year  has  been  spent  in  three  primary  tasks:  i)  ongoing  analyses  of  data  collected  and 
preparation  of  publications  from  the  these  data,  which  are  primarily  from  two  sources  -  cruises  in  2008 
in  the  Bahamas  and  several  cruises  off  Cape  Hatteras  concurrent  with  tagging  of  short-fin  pilot  whales 
(< Globicephala  macrorhynchus)  near  the  Cape  Hatteras  shelf  break. 

The  technical  objectives  were:  i)  to  develop  statistical  models  for  echoes  generated  by  mutli-frequency 
scientific  echosounders;  ii)  to  compare  beaked  whale  echolocation  activity  during  periods  when  our 
echosounders  were  off  vs  on;  and  iii)  collect  data  on  prey  fields  and  physical  oceanography  concurrent 
with  multi-sensor  tagging  of  short  finned  pilot  whales  near  the  Cape  Hatteras  shelf  break;  and  iv)  to 
develop  Bayesian  state-space  animal  movement  models  that  can  be  used  in  the  overall  habitat  models. 

DTAG  Movement  analysis 

In  the  spring  201 1,  we  began  to  investigate  new  modeling  approaches  to  the  dead  reckoning, 
georeferencing,  and  path  estimation  of  DTAG  deployments.  More  accurate  path  estimation  and  explicit 
accounting  of  the  error  inherent  dead  reckoning  will  enable  the  direct  link  between  predator  behavior 
and  concurrent  prey  densities.  Using  deep  diving  pilot  whales  tagged  off  the  coast  of  the  Island  of 
Hawaii,  research  began  applying  Bayesian  state-space  models  to  the  estimation  of  speed  parameters 
and  the  accounting  of  the  error  budget  in  the  dead  reckoning  process.  This  work  continued  into 
summer  2011,  when  the  work  was  conducted  in  collaboration  with  WHOI  faculty  under  the 
supervision  of  Peter  Tyack.  This  work  is  being  conducted  under  the  sub-award  to  Dr.  Pat  Halpin, 
whose  lab  is  responsible  for  assembling  the  habitat  models  that  are  part  of  this  project. 

Statistics  from  Echoes  generated  by  Scientific  Echsounders 

Over  the  last  two  years,  a  Ph.D.  student  at  WHOI  (Wu-Jung  Lee)  has  been  investigating  a  statistical 
approach  to  understanding  the  echoes  generated  by  the  scientific  echosounders.  This  portion  of  the 
project  seeks  to  understand  more  of  what  is  in  the  water  column  from  the  echoes,  i.e.,  adding  to  the 
nonnal  ‘frequency  differencing’  techniques  currently  in  use  to  differentiate  targets. 

APPROACH 

During  the  past  year,  we  have  worked  in  analytical  as  well  as  data  collection  modes.  Our  analyses  and 
synthesis  of  the  prey  field  and  echolocation  data  from  TOTO  has  resulted  a  publication  in  PLoS  ONE 
(Hazen  et  ah,  2011).  St.  Laurent  and  his  group  have  analyzed  much  of  the  physical  data  from  the 
TOTO  experiment,  and  their  results  are  represented  the  M.Sc.  thesis  completed  by  Jay  Hooper  (Hooper 
201 1  -  attached).  Finally,  we  recently  collected  EK60,  ADCP  and  CTD  data  in  the  waters  of  the  shelf 
break  near  Cape  Hatteras,  concurrent  with  attachment  of  DTAGs. 

The  EK60  system  makes  measurements  of  the  returned  energy  from  one  or  more  targets,  which  are 
required  to  obtain  accurate  biomass  estimates  from  the  total  returned  echo  energy  (Foote,  1980a;  Foote 
and  Traynor,  1988).  The  amount  of  backscattered  energy  from  a  single  fish  is  the  backscattering  cross 
section  or  echo  intensity.  If  echo  intensity  is  measured  on  a  log  scale  it  is  called  target  strength  (TS). 
Target  strengths  are  often  measured  during  surveys  (in  situ)  or  predicted  from  trawled  fish  lengths 
using  a  TS-fish  length  regression.  Target  strength  regression  equations  allow  prediction  of  TS  but 
require  large  sample  sizes,  measurements  of  fish  length,  and  often  only  include  one  variable  explicitly. 
For  example,  if  the  tilt  distribution  of  a  school  of  fish  differs  from  the  tilt  distribution  of  the  fish  used 
to  derive  target  strengths,  a  consistent  difference  in  tilt  angle  could  bias  abundance  estimates. 
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Converting  an  acoustic  size  to  a  numerical  size  or  the  total  returned  energy  to  an  acoustic  abundance 
estimate  relies  upon  appropriate  target  strength  values  for  the  population. 

When  fish  echoes  are  too  dense  to  be  counted,  target  strengths  are  required  to  convert  reflected  echo 
energy  to  a  numerical  estimate.  The  linearity  principle  as  defined  by  Foote  (1983)  states  that  the  total 
returned  energy  or  integrated  echo  can  be  divided  by  a  representative  backscattering  cross  section  to 
estimate  fish  abundance.  Fish  lengths  (L)  are  used  in  size-dependent  target  strength  equations:  TS  =  (31 
log  L  +  (30  (1)  where  [31  and  [30  are  parameters  that  vary  among  species  (Love,  1971;  Foote,  1980a; 
Midttun,  1984).  This  target  strength  equation  (1)  explicitly  includes  length. 

Echo  Statistics  in  Mixed  Assemblages 

Also  as  part  of  this  project,  a  PhD  student  at  WHOI,  Wu-Jung  Lee,  has  developed  a  general  theory  for 
the  echo  statistics  due  to  mixed  assemblages  of  scatterers  ensonified  by  scientific  echosounders.  For 
example,  if  there  two  sizes  of  fish  interspersed  and  the  numerical  densities  of  the  two  size  classes  of 
scatterers  are  different,  then  the  echo  statistics  is  distinctly  different  than  the  traditional  "expected" 
Rayleigh  probability  density  function.  Lee’s  general  theory  can  handle  an  arbitrary  number  of  size 
classes  and  numerical  densities;  she  has  drafted  a  paper  on  the  topic. 

For  the  echo  statistics  work,  Lee  has  developed  a  general  theory  for  the  echo  statistics  due  to  mixed 
assemblages  of  scatterers.  For  example,  if  there  two  sizes  of  fish  interspersed  (as  observed  in  a  2008 
experiment  at  Georges  Bank)  and  the  numerical  densities  of  the  two  size  classes  of  scatterers  are 
different,  then  the  echo  statistics  is  distinctly  different  than  the  traditional  "expected"  Rayleigh 
probability  density  function  (PDF).  Her  general  theory  can  handle  an  arbitrary  number  of  size  classes 
and  numerical  densities. 

WORK  COMPLETED 

We  have  analyzed  the  prey  field  data  and  compared  it  with  the  click  activity  recorded  concurrently  by 
the  AUTEC  array  (M3R);  these  analyses  are  contained  in  the  published  manuscript  (Hazen  et  al, 

2011).  The  physical  data  are  being  prepared  for  publication. 

During  the  summer  of  2010  for  the  second  major  portion  of  this  award,  we  worked  as  part  of  an 
interdisciplinary  research  team  to  tag  short  finned  pilot  whales  and  measure  physical  and  biological 
oceanographic  data  at  the  shelf  break  near  Cape  Hatteras,  NC.  We  tagged  8  animals  with  DTAGs  and 
during  these  deployments  we  measured  backscatter  with  38  and  120  kHz  EK60  echosounders,  the 
current  structure  with  a  75  kHz  ADCP,  and  water  mass  characteristics  with  a  CTD. 

In  the  spring  201 1,  we  began  to  investigate  new  modeling  approaches  to  the  dead  reckoning, 
georeferencing,  and  path  estimation  of  DTAG  deployments,  specifically  to  explore  these  movements  in 
areas  where  we  have  prey  field  data.  More  accurate  path  estimation  and  explicit  accounting  of  the  error 
inherent  dead  reckoning  will  enable  the  direct  link  between  predator  behavior  and  concurrent  prey 
densities.  Using  deep  diving  pilot  whales  tagged  off  the  coast  of  the  Big  Island  of  Hawaii,  research 
began  applying  Bayesian  state-space  models  to  the  estimation  of  speed  parameters  and  the  accounting 
of  the  error  budget  in  the  dead  reckoning  process.  This  work  continued  into  summer  2011,  when  the 
work  was  conducted  in  collaboration  with  WHOI  faculty  under  the  supervision  of  Peter  Tyack. 
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RESULTS 


Our  initial  processing  of  the  short-finned  pilot  whale  shows  them  displaying  some  feeding  behavior  at 
depths  where  we  detected  layers  of  scatterers  (Figure  1);  the  whales  also  dove  to  depths  that  did  not 
contain  dense  patches  or  layers  of  prey.  More  analyses  of  these  data  will  follow  soon  with  publication 
expected  within  the  year. 

The  state-space  modeling  work  has  just  begun,  so  we  have  no  results  to  report.  This  work  will  be 
integrated  into  the  final  habitat  modeling,  which  will  be  completed  in  the  final  year  of  the  project. 

For  the  echo  statistics  work,  Lee  has  applied  the  theory  to  broadband  acoustic  data  collected  during  a 
September  2008  cruise  over  Georges  Bank  that  was  not  part  of  the  field  effort  for  the  current  project. 
Both  the  net  samples  and  resonances  in  our  1-6  kHz  system  show  consistent  presence  of  two-size-class 
assemblages  of  fish.  She  has  applied  her  theory  to  the  frequencies  more  relevant  to  the  marine 
mammal  problem,  30-70  kHz.  This  is  driven  by  the  need  to  look  for  cues  in  the  acoustic  echoes  (such 
as  statistical  variability)  for  classifying  in  terms  of  the  sizes  and  types  ("type"  =  zooplankton  vs  squid 
vs  fish,  etc.,  not  at  the  species  level)  that  echo-locating  marine  mammals  might  use  to  determine 
prey.  Her  initial  application  of  the  theory  to  the  data  is  quite  promising  as  can  be  seen  in  Figure 
2.  However,  for  a  rigorous  comparison,  she  needs  to  extend  the  current  theory,  which  is  specific  to 
narrowband  systems,  for  use  with  broadband  data  (specifically,  the  frequency  dependence  of  the 
beampattern  needs  to  be  accounted  for).  She  is  currently  busy  on  that  extension  at  this  point. 

IMPACT/APPLICATIONS 

The  scientific  impact  of  our  accomplishments  are  likely  to  be  significant  because  as  far  as  we  are 
aware,  no  field  program  has  ever  combined  intensive  sampling  of  biological  and  physical  data  in  the 
known  prey  fields  of  a  marine  mammal,  or  any  marine  megavertebrate  for  that  matter.  Furthermore, 
we  are  combining  animal  movement  information  with  these  physical  and  biological  maps  using 
emerging  statistical  behavioral  ecological  models.  The  results  reported  clearly  demonstrate  some 
relationship  between  physical  and  biological  parameters  in  areas  where  beaked  whales  prefer  to  forage 
and  those  they  do  not  frequent. 

RELATED  PROJECTS 

We  have  worked  with  M3R  to  share  the  click  data  from  the  areas  where  we  were  working,  and, 
importantly,  full  audio  data  from  the  time  period  surrounding  our  work  on  the  range  to  assess  any 
potential  effects  our  echosounders  may  have  had  on  the  whales.  A  Master’s  student  at  Duke  has 
analyzed  these  data  for  her  degree,  which  was  completed  in  May  2011  (Vires  2011  and  attached),  but 
she  has  not  yet  completed  publication  of  these  results;  we  will  submit  these  results  for  publication  in 
the  coming  year. 
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Pilot  whale  dive  record  and  sampling  prey  distribution 


Gm_io_2o8a 
July  27,  2010 


250 


500 


14:55 


71-4  ;  H  ■  v-j 


•’  J  •  V*  x  v.*> 


Figure  1.  Time-depth  records  (left)  for  a  diving  short-finned  pilot  whale  carrying  a  DTAG,  and 
water  column  backscatter  (right)  generated  by  the  38  kHz  EK60  system  aboard  the  research  vessel 
following  the  tagged  whale  at  ~ 500  m.  Note  the  dense  layer  of  prey  located  at  the  depth  to  which  the 
whale  is  diving;  also  interesting  is  the  large  apparently  single  target  at  just  less  than  250  m  depth. 
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Figure  2.  Echogram  (right)  and  echo  statistics  (left)  for  assessment  of  target  assemblage  measured 
with  30-70  kHz  echosounders  on  Georges  Bank  during  a  2008  WHOI  cruise. 
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